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ABSTRACT 

Although T Tauri is one of the most studied young objects in astronomy, the nature of its circumstel- 
lar environment remains elusive due, in part, to the small angular separation of its three components 
(North-South and South a-b are separated by 0.68" and 0.12" respectively). Taking advantage of 
incredibly stable, high Strehl, PSFs obtained with Mid-IR adaptive optics at the 6.5 meter MMT, 
we are able to resolve the system on and off the 10 /im silicate dust feature (8.7/im, 10.55/xm, and 
11.86/im; 10% bandwidth), and broad N. At these wavelengths, South a-b are separated by only 
~ 0.3A/D. This paper describes a robust Markov Chain Monte Carlo technique to separate all three 
components astrometrically and photometrically, for the first time, in the mid-IR. Our results show 
that the silicate feature previously observed in the unresolved T Tau South binary is dominated by 
T Tau Sa's absorption, while Sb does not appear to have a significant feature. This suggests that a 
large circumbinary disk around Sa-Sb is not likely the primary source of cool dust in our line-of-sight, 
and that T Tau Sa is enshrouded by a nearly edge-on circumstellar disk. Surprisingly, T Tau Sb does 
not appear to have a similarly oriented disk. 
Subject headings: 



1. INTRODUCTION 

Despite its status as the prototype for young stars, 
the current perception of T Tauri is that it is an ex- 
tremely enigmatic system, perhaps abnormally so. Al- 
though T Tauri w as originally classified as a single star, 
iDvck et 'all (|1982fl discovered that it has an infrared com- 
panion, which ha s neve r been detected in the optical 
(|Stapelfeldt et all (119981) plac e an upper flux limit of 
V ~ 19.6). iGhez et all (|1991h completed an exhaustive 
speckle- image/slit-scan photometric study to construct 
SEDs of the strange infrared companion (hereafter T Tau 
S) along with the original T Tauri (hereafter T Tau N). 
The results showed silicate emission in T Tau N and ab- 
sorption in T Tau S. The incredibly high in frared lumi- 
nosity of T Tau S was enough to convince IGhez et all 
(1991) that it contained its own compact source, and 
flares in their data indicated the presence of an accretion 
d isk around T Tau S . 

iHerbst et all (|1996l) and a followup bv lSolf and Bohm 
(|1999h found perpendicular, molecular outflows: an East- 
West jet is at 23 degrees inclination from the line-of- 
sight and centered on T Tau N, while a North-South jet 
is at 79 degr e es inc lination and centered on T Tau S. 
iKasper et al\ (|2002t ) used integral-field spectroscopy to 
show that Brackett series emission could be constrained 
to a small region around T Tau S, which the authors 
theorized could be indicative of a small edge-on accretion 
disk. 

Compounding the mystery of the infrared companion, 
iRoddier et al\ (|1999l) found t hat T Tau S was non point - 
like and then iKoreskol (|2000l ) resolved T Tau S to be a 

1 The observations reported here were obtained at the MMT Ob- 
servatory, a facility operated jointly by the Smithsonian Institution 
and the University of Arizona. Public Access time is available at 
the MMT Observatory through an agreement with the National 
Science Foundation. 



0.05" binary (here after T Tau Sa and T Tau Sb). When 
observed again bv IDuchene et al\ (|2002ft . T Tau Sb had 
moved si gnificantly, implying th at Sa is a relatively mas- 
sive star IDuchene et ali (|2002f ) were also able to resolve 
Sa from Sb in near-IR spectra, and determined that Sb 
is a pre-main-sequence, early-type M-star with heavy ex- 
tinction and active accretion, while Sa's spectrum is gen- 
e rally feature l ess. 

iBeck et al\ (|2004D presented a study of T Tau's IR 
photometric and spectroscopic variability. The authors 
found that T Tau N is not noticeably variable at K or 
L', while T Tau S varies in the same bands on week-long 
time-scales. While changing accretion rates are usually 
the dominant variability source in classical T Tauri stars, 
the color variability in T Tau S (K-L' vs K) exhibits a 
"redder when faint" phenomenon, which the authors be- 
lieve is best explained by variable extinction. Accretion 
may also be present, but it cannot explain the color vari- 
a bility by itself. 

IDuchene et al\ (|2006f ) were able to use a long baseline 
of observations to astrometrically determine the masses 
and orbital properties of the Southern binary, and the 
results confirmed that T Tau Sa is actually the most 
massive object in the system (N, Sa and Sb have masses 
of -2, 2.73±0.31 and 0.61±0.17M Q respectively). With 
some irony, the star that was originally the prototype 
for the T Tauri star classification, orbits a more massive 
Herbig Ae star. 

Even with the plethora of observations of T Tauri over 
the last 20 years, the nature of the Southern binary is 
still a mystery. Much of this owes to the small separa- 
tions of the 3 components (N-S and Sa-Sb are separated 
by 0.68" and 0. 12" respectively , using the orbital param- 
eterization from lDuchene et all ()2006l )). The fact that Sa 
and Sb are invisible in the optical, and that they can only 
be split with a powerful AO system limits the tempo- 
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ral and spectral ranges over which the whole system has 
been studied. This paper's results, which extend T Tau's 
resolv ed photometric range from 4.7/xm ()Duchene et all 
2003) to 11.86/mi, will improve our ability to study a re- 
solved SED of the Southern binary. Moreover, the spatial 
resolution of the silicate feature definitively establishes 
the stellar source of past unresolved silicate absorption 
detections. 

2. OBSERVATIONS AND REDUCTIONS 

We observed T Tauri in 4 filters (8.7/im, 10.55/im, 
11.86/im, and N-band) with the 6.5m MMT using 
the Mid-IR Array Camera, Gen. 4 (MIRAC4), an 
AO-optimised camera used w ith the MMT ada ptive 
optics system (M MTAO— see iLlovd Hard (|2000( ) and 
iBrusa et all (pOOl )) on Nov. 3, 2006 UT. MIRAC4 is 
a super-sampled, 8-25 /im, 256x256 Si:As array with two 
optical magnifications. We used the high-magnification 
mode (0.055" /pixel) and took 6, 8, 8 and 4 two-minute, 
chop/nod exposures at 8.7 /im, 10.55/im, 11.86/im(10% 
filters), and N-band respectively (Table [T]). We repeated 
these observations for our PSF-star, Beta Gemini (Pol- 
lux). 

At the time of the observations, M1RAC4 was in a 
commissioning, first light run, and had significant elec- 
tronic artifacts. Most of the problems are corrected in 
post-processing using a code written by M. Marengo to 
remove slowly varying channel biases, cross-talk, echos, 
and banding (private communication). Each image is 
then inspected for residual detector effects (weak pat- 
tern noise), and problematic frames are removed (Table 
[T]). Residual detector artifacts are suppressed by median 
combining the images. 

For mid-IR chop/nod observations at the MMT, the 
instrument rotator is left off to stabilize background sub- 
traction by always imaging the same warm, reflective 
surfaces. So to combine our data, we de-rotate all of our 
images by the parallactic angle with cubic spline interpo- 
lation and cross-correlat e the images to ali gn them (stan- 
dard AO reduction as in lClose et all (|2003f 0. The images 
are scaled to the maximum image flux of a centered 30x30 
pixel range to reduce the effects of atmospheric ozone ab- 
sorption, and then median combined (see Figure[T]for an 
example of the reduced 10.55/im image we use). 

The aligned stack of images is used to construct a sigma 
image. In order to robustly calculate the standard devia- 
tion of the image stack, the 3 lowest and 3 highest pixels 
are removed to eliminate outliers (single lowest and high- 
est for broad N). The sigma image is then multiplied by 
a constant to correct for this. 

Because of the spline shifting, and residual detector 
effects, we expect the pixel values and their associated 
errors to be somewhat spatially correlated. However, this 
should only cause a slight smoothing effect that will not 
impact our final results. We parameterize the image into 
two error regions — the star region and the background 
region, and we average the sigma image over the two 
areas to estimate pixel value errors. 

3. ANALYSIS 

3.1. Markov Chain PSF Fitting 

In general, it is possible to super-resolve images at an- 
gular scales far smaller than ~ X/ D given a high enough 



S/N and PSF-stability (jBiller et aLl[2005l ). A binary- 
star deconvolution has 6 free parameters: x%, y\, X2, 
2/2, magi) and ma</2- Adding a third star adds 3 more 
free parameters. Although, mid-IR AO images have ex- 
tremely high and stable Strehls, we find that T Tauri is 
just a little too faint at V-band (MMTAO Strehl starts to 
drop off slightly for V > 11 at 10 /im) to have the near- 
perfec t Strehls our brig hter PSF stars achieve with MM- 
TAO (|Close et aLll2003h . Convolving our PSF stars with 
a 3-parameter (major-axis, minor-axis, angle) Gaussian 
ellipsoid is enough to offset this small effect and allow us 
to work in a mostly photon noise dominated regime. For 
the triple system T Tauri, this gives us 6 position pa- 
rameters, 3 flux parameters, and the 3 PSF modification 
parameters to fit. 

The conventional algorithm to solve a 12-dimensional 
parameter estimation probl em is to do a L evenberg- 
Marquardt x 2 minimization (|Press e t al. 1992). Param- 
eter errors may be estimated from the covariance matrix. 
However, these error bars are unlikely to be even close to 
accurate given the complexity of our problem's parame- 
ter space. Constructing a grid of models and calculating 
relative likelihoods is a much more robust method, but 
a 12-dimensional parameter space is, in this case, too 
large to sample with a grid in a reasonable amount of 
computer time. 

Markov Chain Monte Carlo is the widely accepted 
technique to ci rcumvent the computational unfeasibility 
of mo del grids (jvan Dvldl200llGelman et aZ.ll2003l : iKellvl 
2007). Instead of spending significant CPU time calcu- 
lating model likelihoods in improbable regions, Markov 
Chains follow a random walk where they spend most of 
their time in the most likely areas of the parameter space. 
The direction of the walk is dictated by the Metropolis- 
Hastings algorithm, which allows the chain to spend the 
proper amount of time in each region of parameter space. 
Ideally, the Markov chain has converged when the aggre- 
gate results no longer have any correlation to the chain's 
starting position, and the chain has had time to explore 
every probable part of parameter space multiple times. 
Typically, the beginning of the chain is discarded as a 
"burn-in" phase. To be conservative, we always discard 
the first half of the chain. 

Before implementing the Markov chain, we do a 
Levenberg-Marquardt x 2 minimization to estimate a rea- 
sonable starting value for the Markov Chain. We also 
calculate the covariance matrix from the best-fit, and 
use that to calculate an initial set of jumping conditions. 

For each iteration of the Markov chain, we use a covari- 
ance matrix to draw a proposed jump from a multinomial 
Gaussian distribution. Initially, the covariance matrix 
is a product of the Levenberg-Marquardt best-fit algo- 
rithm, but as the chain progresses, we estimate the co- 
variance matrix from the chain itself. For each proposed 
jump, the model's relative likelihood is calculated and 
compared to the likelihood of the chain's previous step 
(for Gaussian errors, the relative (unnormalized) likeli- 
hood of a model is given by e~ x I" 1 ). If the proposed 
jump is more probable than the previous step, the jump 
is accepted. If the proposed jump is less probable than 
the previous step, it is accepted a certain percentage of 
the time, equal to the ratio of likelihoods of the proposed 
jump and the previous step. 
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When the chain has completed enough iterations, the 
values of each model parameter are binned into a his- 
togram, and fit by a Gaussian. We show these histograms 
for the 10.55^m photometric results in Figure [2] The 
smoothness of the histogram for each parameter is a good 
indication that the Markov Chain has converged. 

At this time, additional quantities of interest may be 
calculated from the Markov Chain. For example, we 
calculate the angular separation of T Tau Sa-Sb us- 
ing the x-y positions of the components, and use the 
known 0.68" angular separation of T Tau N and T Tau 
S (center-of-mass of Sa-Sb, as calcul a ted w ith the orbital 
parameterization of IDuchene et al\ ()2006l )) to measure 
the platescale. 

Using the Markov Chain Monte Carlo technique for 
our PSF-fitting problem has provided us with a sophis- 
ticated and extremely robust tool to estimate the errors 
on our relative photometry measurements. We are also 
able to quantify the degree to which allowing our PSF to 
vary slightly is constrained by the data. The technique 
is general to almost any deconvolution model, however, 
it works best with stable, oversampled PSFs, with rela- 
tively few free parameters. 

3.2. Convergence of the Chains 

We verify that the Markov Chains have converged by 
visually inspecting the histogram of each parameter and 
confirming that the plots are unimodal and effectively 
Gaussian (Figure [2]). We also run our Markov Chain 
Monte Carlo program multiple times with different start- 
ing points and check that our results are unaffected. A 
reasonable reduced x 2 (ie ~ 1) for our best-fit model is 
a good indication that our models and error estimations 
are accurate. Table [2] contains the reduced x 2 01 our 
best-fit models as well as the photometric and astromet- 
ric results for each filter. In Figure [H we compare the 
residuals resulting from a two-star (T Tau N; T Tau S) 
best-fit subtraction, and a three-star (T Tau N; T Tau 
Sa; T Tau Sb) best-fit subtraction. The three-star fit 
has negligible spatial correlations and noise amplitudes 
consistent with photon noise. 

In order to check for systematic effects of our algo- 
rithm, we run blind-recovery tests by constructing fake 
data sets from our resultant parameters and images of 
Beta Gem (our PSF). We add noise to the fake data sets 
consistent with the noise in our real data. The Markov 
Chain technique is always able to return the correct pa- 
rameters within the modeled uncertainties, which conclu- 
sively demonstrates that any systematic effects are lim- 
ited to our slight residual PSF mismatch, and not to any 
algorithmic bias. 

3.3. Photometry and Astrometry 

We do an absolute photometric calibration based 
on known fluxes of ou r PSF star, Beta Gem 
(Hoff mann and Horal Il999). Because of < 10% vary- 
ing atmospheric background/absorption during our ex- 
posures, we assume an absolute photometric calibration 
uncertainty of 10% for each filter, which is a correlated 
quantity between T Tau N, Sa and Sb. The Markov 
Chain technique returns relative errors for each compo- 
nent of T Tau. T Tau N's relative photometric errors are 
dwarfed by the absolute calibration uncertainty, while T 



Tau Sa and Sb's relative errors have an appreciable con- 
tribution to their overall photometric uncertainty. The 
relative errors on T Tau Sa and Sb are almost com- 
pletely anticorrelated (Figure [2ji) , in the sense that if 
Sa is brighter, Sb is fainter and vice versa. 

Astrometric quantities should be invariant across the 
different bandpasses. We use the orbital parameteriza- 
tion from lDuchene et all (|2006l ) along with our measured 
separation of T Tau N and T Tau S (center-of-mass 2 ) to 
determine the plate scale for each observation. This al- 
lows us to measure an angular separation of Sa-Sb, which 
should be constant in each filter. We find that our 10.55 
/um, 11.86 /ira and N-band calculations are consistent 
at the 0.01" level (0.112" ± 0.003", 0.104" ± 0.004" and 
0.121" ±0.002" respectively), and that our 8.7 /xm mea- 
surement is significantly off (0.142" ± 0.004"). However, 
poorer seeing conditions during the 8.7 /im observations 
lead us to believe the longer wavelength results. Varia- 
tions in seeing conditions/Strehl can propagate into the 
strength of the first Airy ring, and because T Tau Sa- 
Sb is aligned with T Tau N's first Airy ring (see Figure 
[Tk andQJ:), separation is a parameter that could be af- 
fected by a small, but noticeable, amount. Systematics in 
our N-band data (see description below) also lend some 
doubt to the accuracy of the broad-band results, which 
could explain why ou r error bars ap p ear to be slightly 
underestimated. The IDuchene et al\ ()2006l ) orbital pa- 
rameters predict a separation of 0.119" ± 0.004" at the 
time of our observations, which is consistent with our 
measurements at the 0.01" level. 

The position angle of the Southern binary is well con- 
strained and consistent across the filter set (307.4 ±1.1. 
306.5 ± 1.6, 308.4 ± 1.1 and 311.1 ± 0.9 degrees at 
8.7/n m, 10.55jnm, 11.86/xr n and N-band respectively). 
The IDuchene et al\ ([20061 ) orbital parameters predict a 
position angle of 306 ± 13 degrees at the time of our ob- 
servations, which is consistent with our measurements, 
although our data place a much tighter constraint. 

In Figure |4] our astrometry resu lts are compare d to th e 
orbital motion measurements of Scha efer et all ([2006) . 
Our results are generally consistent with the predicted 
orbits. While there is some preference for orbits shorter 
than 40 years, our results are meant more as a demon- 
stration of fitting accuracy than as an orbital constraint. 

Our N-band results are somewhat questionable because 
we only had data from two different chop-nod sets (and 
only six usable images total). This means we could not 
adequately suppress the detector's correlated noise by 
median combining the images, and as a consequence, sys- 
tematics may have dominated our results at N-band more 
than for the other filters. Our best-fit reduced \ 2 a t N- 
band is 4.43 (as opposed to 1.24, 0.82 and 1.37 at 8.7/xm, 
10.55/im and 11.86/im respectively). 

4. DISCUSSION 

Figure [5] shows the fluxes and overall flux errors of 
T Tau N, Sa and Sb in the three narrowband filters. 
These results show that t he silicate absorpti on in T Tau 
S, originally observed bv iGhez et al\ (|1991h . originates 
entirely from T Tau Sa. T Tau Sb's silicate SED has 

2 The center-of-mass is calculated by taking our measured sepa- 
ration of T Tau S a-Sb and assuming a mass ratio for T Tau Sa-Sb 
as determined bv IDuchene et al\ (120061 ') ■ 
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large enough photometric uncertainty that we cannot say 
whether it has a small emission or absorption feature. 
However, we can say with certainty, that it does not have 
as dramatic a silicate feature as T Tau Sa. T Tau N also 
h as a negligible silicat e feature. 3 

I Walter et aH (|2003| ) used STIS spectra to infer the 
presence of circumbinary structure obscuring both ob- 
jects of the T Tau S binary. However, a circumbinary 
structure cannot be the major source of silicate absorp- 
tion since the absorption is only towards Sa. Our obser- 
vations indicate the presence of an edge-on circumstellar 
disk around T Tau Sa. The dense, optically thick mid- 
plane of a disk can obscure the star (causing the high ex- 
tinction necessary to completely hide a ~ 2.7 solar mass 
star shortward of H-band) and create the silicate absorp- 
tion we observe in an otherwise featureless spectrum. T 
Tau Sb is less red, less obscured a nd has less (if any) sil - 
icate absorption than T Tau Sa. IDuchene et all (2002) 
spectra indicate that T Tau Sb is an MO star with heavy 
extinction and active accretion, implying the presence of 
a disk. Combined with our null detection of silicate ab- 
sorption, it is likely that T Tau Sb has a non-edge-on 
disk. 

In general, circumstellar disks in tight binaries are 
tidally alig ned on shor t time -scales, but iJensen et all 
(2004) and lMonin et al\ ()2006l ) have found that systems 
with three or more stars tend to have misaligned disks. 
These previous works used polarimetry, and could only 
separate binaries with sep > 100AU. By resolving sili- 
cate features with 10 /im AO and super-resolution tech- 
niques, we can push this limit to ~ 15 AU where tidal 
forces are stronger and disks are truncated. 

5. CONCLUSIONS 

Markov chain super-resolution is a useful technique for 
taking full advantage of highly stable, diffraction limited 
images from ground-based AO systems or space-based 
telescopes. The MMTAO system, along with MIRAC4's 
supersampled detector, is uniquely capable of producing 
these images in the mid-IR at a 6.5-meter class telescope. 

We have used MMTAO to image the famous T Tauri 
triple system and were able to split the 0.11" Southern 
binary on and off the silicate feature (8.7/^m, 10.55^m, 



3 IGhez et ali ((T99lh were able to see an emission feature from T 
Tau N, but their SED shows that the peak of the emission is proba- 
bly fairly localized to 9.7/mi, which would explain our null result at 
10.55/mi. The amorphous silicate feature generally peaks at 9.7 (im, 
however atmospheric ozone absorption makes this wavelength hard 
to observe from the ground. Strong features and crystalline sili- 
cate can still be detected in the 10.55/^m filter. Our N-band results 
(using N — [8.7/^m-t- 10.55^m + 11.86/im] as a proxy for a "9.7/im" 



11.86/im and broad N-band). While previous unresolved 
data from the Southern binary show strong silicate ab- 
sorption (explaining its high extinction), we have deter- 
mined the source of the absorption is entirely in front 
of T Tau Sa, a ~2.7 solar mass Herbig Ae star with an 
otherwise featureless spectrum. 

Our results indicate the presence of an edge-on circum- 
stellar disk around Sa, corroborating previous theories 
resulting from jet orientation, differential extinction of Sa 
and Sb, narrowly constrained Brac kett series em i ssion , 
and warm, narrow CO absorption dHerbst et al\ ["1 996; 
Solf and Bohdll999t iKasper et a/.ll2002l : IDuchene et ali 



2005T) . 



T Tau N has been observed to have silicate emis- 
sion (IGhez et all |l991) and has a perpendicu l ar jet 
orientation to the Southern jet (He rbst et ali 119961 : 
ISolf and Bohmlir999D . which indicates the presence of a 
face-on disk. Sb has been modeled with an accretion disk 
and the lack of silicate absorption in our results means 
that it also has a non-edge-on disk. The fact that disks 
are misaligned in such a tight (P ~ 20 year for Sa-Sb) 
system is surprising considering that tidal forces should 
align them on short timescales. As usual, the prototype 
for young stars has provided another surprise for star 
formation models. It remains to be seen whether other 
tight binaries/triple systems have similar anomalies. 
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filter) corroborate the fact that emission/absorption are peaked at 
9.7/im (although, as was mentioned in the previous section, the 
accuracy of the N-band results is questionable). Sb has a higher 
flux at N-band than any of the narrow-band filters, indicating that 
Sb could have significant emission at 9.7 /im. North also has a rea- 
sonably strong N-band point, which suggests that it also probably 
has silicate emission. 
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TABLE 1 

Observations of T Tau (Nov. 3, 2006) 



parameter 


8.7 fim 


10.55/^m 


11.86/im 


N 


Chop-Nod Sets a 


6 


8 


8 


4 


On-Source Time (sec) 


720 


960 


960 


480 


Usable Images 


20 


25 


20 


6 


Usable On-Source Time (sec) 


600 


750 


600 


180 



a Each chop-nod set comprises four images with 8" chops and 6" nods 



TABLE 2 

Photometric and Astrometric Measurements (Nov. 3, 2006) 



parameter 



;.7^tm a 



10.55/im a 



11.86Aim a 



N b 



T Tau N (relative mag) 

T Tau Sa (relative mag) 

T Tau Sb (relative mag) 

T Tau Sb-Sa (delta mag) 

T Tau N (Jy) 

T Tau S (Jy) 

T Tau Sa (Jy) 

T Tau Sb (Jy) 

Sa-Sb Separation (arcsec) 

Sa-Sb PA (degrees) 

best reduced \ 2 



0.000 ± 0.002 
-0.029 ±0.011 
1.977 ±0.075 

2.005 ± 0.087 
6.87 ±0.69 
8.17 ±0.82 

7.06 ±0.71 
1.11 ± 0.13 
0.142 ±0.004 
307.4 ± 1.1 
1.24 



0.000 ± 0.001 
1.134 ±0.058 
1.639 ± 0.092 
0.498 ± 0.151 
7.51 ±0.75 
4.30 ±0.43 
2.64 ±0.30 
1.66 ±0.22 
0.112 ±0.003 
306.5 ± 1.6 
0.82 



0.000 ±0.001 
0.244 ±0.047 
1.215 ±0.117 
0.966 ±0.165 
7.43 ± 0.74 
8.36 ± 0.84 
5.93 ±0.65 
2.43 ±0.36 
0.104 ±0.004 
308.4 ± 1.1 
1.37 



0.000 ±0.002 
0.634 ± 0.030 
0.702 ± 0.033 
0.064 ± 0.063 
7.40 ±0.74 
8.01 ±0.80 
4.13 ±0.41 
3.88 ± 0.39 
0.121 ±0.002 
311.1 ±0.9 
4.43 



a 10% filters 
b 8-13 pm 
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Fig. 1. — For all three graphics, North is up and East is left. The images are 1.7" on a side and have a logio stretch. 
Left: The 10.55/xm median combined image of T Tauri showing T Tau N and T Tau S. Here, T Tau S appears to be unresolved. 
Center: The 10.55/im median combined image of Beta Gemini (the PSF used in this paper). 

Right: A schematic o f the T Tauri system, with 10.55/im photometry and astrometry as derived in this paper (the N-S separation is from 
IDuchene et a l. (2006)). The PSF used in the schematic is scaled down for aesthetics. 
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magnitudes (relative to T Tau N) T Tau Sa (Jy) 



Fig. 2. — Photometric results of the Markov Chains for the 10.55^im filter. The T Tau N, Sa and Sb plots are relative magnitudes scaled 
to the mean result for T Tau N. The fourth plot is a scatter plot of fluxes for T Tau Sa vs. T Tau Sb. The —1 slope demonstrates that 
the T Tau Sa and T Tau Sb flux errors are anticorrelated (if Sa is brighter, Sb is fainter and vice versa). The fluxes in the fourth plot are 
relative, and thus ignore the 10% photometric calibration error assumed in the rest of the paper. 
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Two Star Best-Rt I hree Star Best-Fit 




Fig. 3. — Levenberg-Marquardt best-fit residuals at 10.55/im with two-star(left) and three-star models (right). The three star fit has 
residuals at the photon noise floor. 
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Fig. 4. — Figure modified from Schacfcr et al. (2006) 
Astrometric solutions to the T Tau Sa-Sb binary (Sa is in the center) with this paper's results (2006.9) for f0.55/im (red), 11.86/im (green), 
N-band (blue) overplotted. There is some preference for orbits with a period shorter than 40 years. 




Fig. 5. — The figure shows photometry for 8.7/mi 10.55^tm and 11.86^tm for T Tau N, Sa, and Sb. The curves drawn through the points 
are intended as a visual aid. T Tau Sa has a large absorption feature that is absent from the other stars, which indicates the presence of an 
edge-on protoplanetary disk. Since the other stars lack a similar feature, it is likely that the disks in the T Tauri system are misaligned. 



